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Abstract 

The high magnetic concentration Ising antiferromagnet Feo.93Zno.07F2 does not 
exhibit the severe critical scattering hysteresis at low temperatures observed in all 
lower concentration samples studied. The system therefore provides equilibrium neu- 
tron scattering line shapes suitable for determining random-field Ising model critical 
behavior. 



It is well known that a phase transition occurs [p]] for the d — 3 random-field Ising 
model (RFIM). However, critical behavior investigations using the experimental realizations, 
anisotropic randomly dilute antiferromagnets such as Fe x Zri\- x F2 in applied uniform fields, 
have been very challenging, primarily as a result of severe hysteresis in scattering line shapes 
below the transition T C (H). This is most evident when comparing data taken upon cooling in 
a field (FC) and heating in the field after cooling in zero field (ZFC). It has been shown [0, |3| 
that, at large dilution, the system breaks into weakly interacting domains as T — > T C (H) 
after ZFC to establish long-range order. For x < 0.8, vacancies are so numerous that 
domain walls form with little energy cost. This effect obscures the decrease of the order 
parameter within domains due to thermal fluctuations since the Bragg intensity vanishes 
with domain formation. In the present study, we have overcome this problem by employing 
a magnetically concentrated crystal, Feo.g3Zrio.07F2, in which domains cannot form without 
the energy penalty of breaking many magnetic bonds, as in the original ferromagnetic RFIM. 

The scattering measurements were made at the Oak Ridge National Laboratory High 
Flux Reactor using a two-axis spectrometer configuration. We used the (0 2) reflection of 
pyrolytic graphite (PG) at an energy of 14.7 meV to monochromate the beam. We employed 
three different collimation configurations. The lowest resolution is with 70 min of arc before 
the monochromator, 40 before the sample and 40 after the sample. We also took scans with 
20 min of arc before and after the sample and with 10 min of arc before and after the sample. 
PG filters were used to eliminate higher-order scattering. The carbon thermometry scale was 
calibrated to agree with recent specific heat results for the H = transition. The field 
dependence of the thermometry was also calibrated and T C (H) — T c (0) is consistent with the 
specific heat data. The rounding of the transition occurs only for \t\ < 2 x 10 -3 , allowing 
the critical behavior to be probed reasonably well. 

Four transverse ZFC scans across the antiferromagnetic (10 0) Bragg point with fine 
collimation are shown in Fig. 1. Two scans are at 0.13 and 0.31 K above T C (H) = 70.62X 
and the others are 0.14 and 0.20 K below. The contrast between the scans above and below 
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T C (H) is striking. All scans taken at H = and those below T C (H) for H > can be 
reasonably fit to the mean-field Lorentzian line shape 



= - + M s 2 5(q) , (1) 
q z + k z 

where k(T) is the inverse correlation length for fluctuations and M s is the staggered magne- 
tization. The mean-field RFIM prediction includes a squared Lorentzian term 

s M = ^+(^ + W • < 2 > 

but the H > and T < T C (H) fits allow no significant contribution from this term. Above 
T C (H) Eq. [I] does not fit the data well and Eq. |2] is more appropriate. The fitted behavior 
for the staggered susceptibility, x{T), and k{T), obtained from the Lorentzian term, does 
not behave consistently as explained below. We therefore emphasize that these equations 
are not adequate, perhaps because the RFIM exponent rj, a measure of the deviation from 
mean-field behavior, is quite large JTJ. 

The q = Bragg intensity is very large just below T C (H) and vanishes just above, in 
accordance with Monte Carlo results || that (3 is close to zero. The small hysteresis observed 
in the Bragg intensity upon FC and ZFC is consistent with the extinction effects. 

The \q\ > line shapes also show a dramatic change at T C (H) with no hysteresis. For the 
previous studies at x = 0.52 |2| and 0.46 ||, the \q\ > intensities grow much more rapidly 
as T — > T C (H) and, concomitantly, the Bragg scattering exhibits an anomalous decrease in 
intensity, effects attributed || to domain formation which are apparently absent for x = 0.93. 

k(T) and x{T) are shown in Fig. 2 and Fig. 3, respectively, for both H = and H = 7T. 
In both (T) decreases to approximately the resolution limit of 2 x 10 3 reciprocal 

lattice units. At H = 0, the values of k(T) and xCO ^ well to random-exchange behavior 
H above and below T C (H), with v = 0.71 ± 0.01 and 7 = 1.35 ± 0.01, respectively. For 
H = 7 T, the values v = 0.90 ± 0.01 and 7 = 1.72 ± 0.02 for T > T C (H) are consistent 
with earlier experiments 0, simulations || and theory p. The disconnected staggered 
susceptibility exponent 7 = 3.0 ±0.1 is also in reasonable agreement with previous results 
0. However, the values of k(T) and x(T) for T < T C (H) are inconsistent with any power 
laws and fits are consequently not shown. The results for k(T) and x(T) for T < T C (H) are 
very dependent on the scattering line shapes used in the analysis and further discussion will 
be deferred to a more extensive report. The values for k(T) and x(T) are independent of 
the cooling procedure, implying equilibrium behavior. 
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Figure 1: The logarithm of the neutron scattering intensity vs. q just above and below 
T C (H). 
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Figure 3: The logarithm of the staggered susceptibility \ vs. T for H = and 7 T. 
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